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We propose a microscopic model that describes the mag- 
netic behavior of the mixed-spin quantum systems R2BaNi05 
(R= magnetic rare earth). An evaluation of the properties 
of this model by Quantum Monte Carlo simulations shows 
remarkable good agreement with the experimental data and 
provides new insight into the physics of mixed-spin quantum 
magnets. 

PACS numbers: 75.25. -fz, 75.10.Jm, 75.50.Ee 

Low-dimensional (low-D) quantum magnetism remains 
at the forefront of condensed matter research for al- 
most two decades. This is primarily because this field 
deals with simple models of magnetism that demonstra- 
te a broad spectrum of complex quantum-mechanical 
phenomena of a general and fundamental nature. The 
simplicity of the Hamiltonians involved, as well as the 
low dimensionality, often allow accurate theoretical and 
numerical treatment. On the other hand, the discov- 
ery of real low-D materials allows experimental stud- 
ies for direct comparison with theory on the quantita- 
tive level. A particular area of great recent interest is 
the co-existence of "quantum" and "classical" proper- 
ties in systems composed of weakly coupled quantum 
spin chains. One good example is the interplay between 
continuum and spin-wave dynamics in gapless quasi- 
one-diipjensional (quasi-lD|) Heisenberg systems such as 
KCuFsEI and BaCu2Si2 07l3. At the focus of the present 
work is a more complex phenomenon: the seemingly 
paradoxical co-existence of 3D magnetic long-range order 
and ID quantum gap excitations in rare earth nickelates 
with tie. general formula R2BaNi05 (R= magnetic rare 
earth)Eiy. 

R2BaNi05 species have two types of spin carriers. 
5 = 1 Ni^+ ions form distinct antiferromagnetic chains 
running along the a axis of the crystal structure. If these 
chains were perfectly isolated, they would be a classic ex- 
ample of a Haldane antiferromagnet with no long-range 
order even at T^— and a gap in the magnetic exci- 
tation spectrumO. Long-range ordering that occurs in 
R2BaNi05 at low temperatures is driven by the mag- 
netic s = I R'^"'' ions, positioned in-between the chains. 
3D magnetic order involves both the R"^+ and Ni^+ spins. 
What is particularly interesting though, is that Haldane- 
gap excitations associated with the Ni-chains persist in 
the ordered state and co- exist with conventional order- 
parameter excitations (spin waves). By now, a wealth of 
experimental data|-|(j«-|Several R2BaNi05 compounds has 
been accumulatedQ'oB. Most theoretical work, however, 
was based on a simple Mean Field (MF)q or Random 



Phase Approximation (rpa)I!I treatment of the interac- 
tions between Ni-chain and rare earth subsystems. The 
existing numerical studies were aimed at calculating the 
"bare" properties of the Ni-chains, needed to complete 
the MF/RPA equations^. While the MF/RPA approach 
turned out to be an extremely useful tool in understand- 
ing the basic physics involved, it has numerous intrin- 
sic limitations, particularly for dilute (Ri_2:Y2:)2BaNi05 
systemscl, or in the case of strong Ni-R coupling. In 
the present work we abolish the MF/RPA framework 
and construct a microscopic model for R2BaNi05 com- 
pounds. We then employ this model in a systematic first- 
principle numerical study of these interesting materials. 

In order to construct the microscopic model we will 
use the available experimental results to design an ap- 
propriate description for the magnetic ions and interac- 
tions in the system: (i) A good understanding of the 
Ni-subsystem can be drawn from the known properties 
of Y2BaNi05, a material where the magnetic rare earths 
have been replaced by non-magnetic Y'^"'" . The Yttrium- 
nickelate is an almost perfect physical realization ofj-a 
Haldane gap antiferromagnet. Neutron scattering datatHi 
show a large gap A ^ lOmeV and well characterized Hal- 
dane excitations. This suggests that the intrachain Ni-Ni 
exchange coupling in the chains running along the a-axis 
has to be J '--^ 25meV. Furthermore, there is no evidence 
of 3D long range order even at the lowest temperatures. 
Any residual direct coupling between the chains, both in b 
or c direction, is smaller than the critical value necessary 
to establish AF order at low enough temperaturesEJE3. 
This critical value depends quadratically on the gap and 
it is expected to be large in this family of nickelates due to 
the magnitude of the gap. Therefore we will consider the 
direct Ni-Ni interchain coupling irrelevant for a physical 
description of the material. 

(ii) When Y, located between the chains, is completely 
substituted by a magnetic rare earth R, the system orders 
antiferromagnetically. The Neel temperature is typically 
smaller than the gap energy. For example, in the_.case 
R=Nd the Tn = 48K and A - llmeV (- 127K)EI. In 
the paramagnetic phase (T > Tn), the Ni-chain excita- 
tions in R2BaNi05 are virtually indistinguishable from 
those in the Yttrium compound at the same tempera- 
ture, implying that very little should change in the model 
of the S ~ 1 chains while the ordering should be driven 
by some kind of exchange Ni-R that induces an indirect 
coupling between the chains. 

(iii) In the magnetically ordered phase, the staggered 
magnetization of the Ni sublattice has as T ^ satura- 
tion values of 1.0 — 1.6fj,B per ion in all the R2BaNi05 
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studied to date^M'H These values are clearly smaller 
than the classical result 2hb- This observation sig- 
nals that the model describing these compounds has to 
properly retain the quantum fluctuations in the ordered 
phase. 

(iv) Unlike Ni-R interactions, that are vital in order to 
induce a static magnetization on the Ni sites below Tn, 
direct R-R magnetic coupling can be disregarded in our 
model. Indeed, the ordering temperatures in R2BaNi05 
compounds are typically several tens of Kelvin. Direct 
coupling strength between rare earths in insulators is 
much weaker, usually of the order of one Kelvin. This 
is due to the fact that in rare earth species the mag- 
netic /-electrons are strongly localized, what prevents 
efficient superexcliaMe couplingtB. Unlike the case of 
Y2_2;Ca2;BaNi05llII"ll^, where doping with Ca introduces 
hole carriers in the Ni-chains, this localized nature of f- 
electrons in the rare-earth keeps the S=l chains free from 
charge carriers. 

(v) In the crystal structure of R2BaNi05, the site- 
symmetry for R'^"'" is very low. As a result, for each 
rare earth ion the degeneracy of its magnetic multiplet 
is lifted by crystal field effects. To investigate the low- 
energy part of the spectrum we only need to consider 
the lowest-energy orbital levels. In the case of the rare 
earth beina-a. Kramers ion with half-integer total angular 
momentumll3, we shall model the magnetic rare earths by 
effective s = ^ (pseudo)spins. In order to decide about 
the type of coupling between the Ni and the R sublattices, 
we rely on the experimental observation of the absence of 
dispersion in the crystal field excitations associated with 
R^+ ionsD. This fact suggests that the coupling between 
the R and Ni sublattices is very anisotropic and can be 
approximated by an Ising-type term in the Hamiltonian. 

With all the ingredients above, we propose now the 
following Hamiltonian: 

^ = Si^2jSi+i,2j + Jc^^Sl2j{sl2j-l + sf,2j + i) 

ij ij 

with J and Jc both positive. The first term in is a 
5 = 1 Heisenberg model along the Ni chains, the second 
is an Ising-like coupling with strength Jc between the 
5 = 1 Ni ions and the s = 5 magnetic moment in the 
rare earth. The index i runs along the chain direction 
and J in the direction perpendicular to the chains. A 
sketch of the model is shown in Fig. (|l|). When Jc = 0, 
H should reproduce the physics of the Y2BaNi05, i.e. 
basically independent S = 1 chains. 

We have studied this model numerically with the 
Quantum Monte Carlo (QMC) methodjising a quantum 
cluster algorithm, the Loop AlgorithmEil. This method 
allows for an efficient sampling of the configuration space 
giving a specific prescription of how to perform global 
updates in quantum systems. Such updates involve clus- 
ters of spins with a size of the correlation length. The 
hamiltonian we studied is not frustrated and does not 
show the sign problem. All the Boltzmann weights can 



be taken positive after the conventional rotation around 
the z-axis of all the spins in one of the two sublattices of 
the S=l system. The difhculty in reaching extremely 
low temperatures is a drawback of the QMC method 
when one is interested in characterizing exotic ground 
states, mostly near quantum critical points or strongly 
frustrated systems. None of these cases affects our calcu- 
lations since our main goal is to compare with experimen- 
tal results at temperatures, as we will see, perfectly ac- 
cesible for QMC. We have measured various magnitudes 
that can be directly compared with the available experi- 
mental data and which, in some cases, are not accesitJg 
for other approaches like the MF/RPA approximationQ-lD 
or the zero temperature Density Matrix Renormalization 
Group (DMRG) calculationstl. These magnitudes are: (i) 
Magnetizations in both sublattices (ii) Neel temperatures 
(iii) Correlation lengths and gaps in both subsystems. 
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FIG. 1. Microscopic model for R2BaNi05. The big circles 
denote S = 1 Ni spins and the small circles denote the s = i R 
spins. J and Jc are the Ni-Ni and the Ni-R exchange couplings 
respectively. 

Before discussing the results, we present a brief tech- 
nical description of the numerical work and data treat- 
ment. QMC simulations are performed in finite lattices 
therefore an exhaustive finite size analysis of the results is 
necessary to compare with real experiments. We used pe- 
riodic boundary conditions. The magnetization (M) was 
computed by directly measuring |M| in order to avoid the 
averaging between configurations around the two ground 
states of the model in the ordered phase, and by using the 
relation M — Ximj^^ao l-^l- We considered for this ex- 
trapolation LxL lattices of size L=8, 16,24,32,48,64 spins. 
The Neel temperature was then determined by using the 
Binder parametern3 which is the fourth cumulant of the 
order parameter distribution. At Tn this cumulant is 
independent of the size of the system, apart from sub- 
dominant corrections. We have observed that for values 
of Jc smaller than ^ O.IJ these corrections are due to the 
anisotropy in the scaling (i.e. different number of spins 
in a and c directions). We have checked this point by 
studying systems with different sizes La — 32,48 spins 
and Lc = 24, 32, 48 spins. For values of Jc ^ 0.3 and 
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the precision needed for comparison with the experiment, 
T/v can be computed by extrapolation from isotropic lat- 
tices. The statistical error bars are always smaller than 
the symbol sizes in all the figures presented. We mea- 
sured all the magnitudes in 10^ Monte Carlo steps after 
thermalization. We have used units = = h = 1. 

The staggered magnetizations obtained from the QMC 
computations for both the Nickel (A/ni) and rare earth 
(Mr) sublattices are in good agreement with the exper- 
imental results for Nd2BaNi05a for a suitable choice of 
the model parameters in the Hamiltonian. For instance, 
for a value of the transverse coupling Jc = 0.31J, we 
have that MNi(T' — > 0) = 0.79 that corresponds to the 
value 1.6/iB observed experimentally. The Neel temper- 
ature for this value of Jc is Tn = 0.163 J. We can take 
the gap in a single S = 1 chain, A = 0.410 J and com- 
pute the ratio between the two main energy scales in 
the system r = ip-. We observe that r is very similar 
in the material (r=2.59) and in the model (r=2.51). In 
Fig. (H) we present the temperature dependence of the 
staggered magnetizations in the S = 1 (open squares) 
and the s — ^ (open circles) sublattices, A^ni and AfR 
respectively, for the selected value of Jc obtained from 
the QMC simulations. The magnetization results have 
been rescaled with respect to the maximum saturation 
values, i.e. 1 for the S = 1 system and 0.5 for the s = | 
system. The results can be compared with the experi- 
mental magnetic moment in the Ni (filled squares) and 
Nd systems (filled circles)a. Note the good agreement 
between the theoretical and experimental results at low 
temperatures. We also observe that the staggered mag- 
netization in the S = 1 sublattice is nearly temperature 
independent for T < At T=40K, theNdaNiBaOs 

undergoes a magnetic structural transitioned and there- 
fore we expect only qualitative agreement with our model 
at temperatures near Tjv. 
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FIG. 2. A/ni (squares) and A/r, (circles) vs. T. Open sym- 
bols are the results obtained with the model Hamiltonian for 
Jc=0.31J and filled symbols are the experimental data. Both 
magnetizations are given in units of /xs and have been rescaled 
with respect to the saturation values 1 and 0.5 respectively. 

Both sublattices behave qualitatively different in the 
ordered phase. While in the Ni sublattice the reduction 



of the staggered magnetization by purely quantum fluc- 
tuations is important, the R sublattice is fully saturated 
to the classical value A^r = 0.5 at T=0. 

In order to understand the origin of the ordered phase 
of this quantum model we have studied the relation be- 
tween the magnetizations in both sublattices for different 
values of Jc- In Fig. (^(a) we show a plot of the QMC 
results at T < for Mr as a function of Mj^n for var- 
ious values of Jc- We observe that the behavior of the 
magnetizations is well reproduced by the function: 



A/r = A/otanh(a/3A.fNi) 



(1) 



where (3 = l/ksT, Mo is the effective moment of the 
rare earth ion and a is a linear function of Jc in the 
range 0.09 < Jc < 0.5J as observed from our QMC re- 
suhs. Eq. (0) can be obtained by a mean-field approach 
as shown by Zhcludcv et ala. These authors considered 
that the behavior of the RaBaNiOs compounds in the 
ordered phase could be described in terms oi a S — 1 
chain of Ni ions in a staggered magnetic field induced by 
the s = i magnetic rare earth ions and, reciprocally, the 
otherwise free spins s = 5 see the mean field produced 
by the neighboring 5 = 1 chains. 
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FIG. 3. QMC results for (a) Mr vs. Mni for various 
values of the couplingconstant Jc- The solid lines correspond 
to the function Eq. (|l|). (b) Mni-ivs. the 'effective' field on 
the Ni-subsystem induced by A/rQ for vaiiious T values. The 
squares correspond to the DMRG resultaS at T=0. 

The excellent agreement between this mean-field ap- 
proach and our QMC results indicates that our micro- 
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scopic model also fulfills the property of viewing the ef- 
fect of the staggered magnetization on the rare earth ions 
as an efi^ective magnetic field on the Ni subsystem and 
viceversa. 

In Fig. (|^)(b) we show the relation between A/ni and 
the staggered magnetic field induced by Mr obtained by 
our QMC calculations for various T values. The extrapo- 
lation T ^ reproduces the DMRG results for one S = 1 
chain in a. staggered magnetic field at T = obtained by 
Yu et aln. Since the R sublattice is fully polarized at 
T = no rescaling is necessary to compare our data with 
the DMRG results. The quality of the fits in Fig. (|)(a) 
and the data collapse in Fig. (||)(b) suggests a technical 
procedure to extrapolate QMC data down to zero tem- 
perature. This procedure is obvious for magnetizations 
but could be extended for other thermodynamic magni- 
tudes for Hamiltonians similar to the one proposed here. 

In the context of the effective field interpretation, an 
analysis of the temperature dependence of Mr obtained 
with our QMC simulations indicates that though Mr sat- 
urates to the classical value at T=0 (see Fig. (||)), as ex- 
pected if we consider this subsystem as a sublattice of 
free spins in the presence of an external field, this behav- 
ior cannot be reproduced by a purely classical mean field 
model in the R sublattice, namely mean field just pro- 
portional to the magnetization in the R sublattice and 
independent of the Ni magnetization. This fact is sig- 
naling that quantum fluctuations in the Ni sublattice are 
essential to describe the magnetization in the R sublat- 
tice. 




FIG. 4. Inverse correlation length vs. temperature for 
the Ni-chain subsystem (circles), the R-subsystem along a 
(squares) and for a single 5* = 1 chain (triangles). Solid lines 
correspond to the law ^"^ ~ (T - Tm)^ ■ At higher T, the 
Ni-subsystem behaves like a Haldane 5* = 1 chain. 

We discuss now the spatial spin-spin correlations in the 
paramagnetic phase. In Fig. (|^) we plot the inverse of the 
correlation length along the Ni chains (circles) and 
along the R chains parallel to the Ni chains (squares) for 
T < Y- For comparison, we present in a single S — 1 
spin chain (triangles). In this range of temperatures it is 
well known that the correlation length in the = 1 chain 



is ^ = ^ where c is the spin-excitation velocity. The 
correlation length is computed in all the cases by fitting 
the z-camponent of the spin-spin correlation function to 
the lawO: 



(2) 



(A and rj are fitting parameters) directly in a lattice of 
64x64 spins and for Jc = 0.18 to keep and A well 
separated in energy. Near T/v, the critical modes become 
gapless and we observe how the gap closes in both sub- 
systems following a law = K{T — Tjy)^ . K is the 
only free parameter in the fit presented in Fig. (^ as a 
solid line, since T/v is computed independently. As the 
temperature increases, the correlation length in the Ni 
subsystem approaches the correlation length of a single 
chain. 

The behavior of 77 is completely different in both sub- 
lattices. While in the R sublattice the correlation func- 
tion (SqSi) can be fitted with = for all the tempera- 
tures showed in Fig. (H), in the Ni sublattice r] approaches 
the value 0.5 in Eq. (0) as the temperature is reduced in 
the paramagnetic phase. We have analized this behavior 
and observed that the temperature dependence ofrjiT) in 
the spin-spin correlation function Eq. (^) is very similar 
to that of 77 in a single chain in which Eq. (^ shows a slow 
cross-over from 77 = i at T=0, more precisely {SqSI) is 
then described by the. modified Bessel function K(B^, to 
77 = when T > Z\c3. The temperature dependence of 
77(T) is shown in Fig. (||). 
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FIG. 5. Temperature dependence of the fitting parameter 
ri{T) in Eq. for the Ni-subsystem . 

The picture arising from the study of the spin-spin 
correlation functions is that the modes becoming gapless 
at Tn can be described at a mean field level. Moreover, 
since both the value of the correlation length and the 
power law corrections to the exponential behavior are 
very similar for the Ni subsytem and for an independent 
5 = 1 chain, we can conclude that the modes remaining 
gapped in the system are very much like the conventional 
Haldane excitations. 

In summary, we have shown that QMC calculations 
based on the proposed microscopic model are able to re- 
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produce the observed behavior of R2BaNi05 mixed-spin 
quantum antifcrromagncts remarkably well. They pro- 
vide a solid numerical basis for the MF/RPA model, and 
make additional predictions regarding the magnetic na- 
ture of both spin species that go beyond the simplistic 
picture of a MF/RPA model. We hope that the proposed 
approach will be particularly useful in the study of di- 
lute (Ra;Yi_j;)2BaNi05 compounds in the limit of small 
X, where the MP model becomes totally inadequate. 
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